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In this paper we study a 3D lattice spin model of CP 1 Schwinger-bosons coupled with dynamical 
compact U(l) gauge bosons. The model contains two parameters; the gauge coupling and the 
hopping parameter of CP 1 bosons. At large (weak) gauge couplings, the model reduces to the 
classical 0(3) (0(4)) spin model with long-range and/or multi-spin interactions. It is also closely 
related to the recently proposed "Ginzburg-Landau" theory for quantum phase transitions of s = 1/2 
quantum spin systems on a 2D square lattice at zero temperature. We numerically study the phase 
structure of the model by calculating specific heat, spin correlations, instanton density, and gauge- 
boson mass. The model has two phases separated by a critical line of second-order phase transition; 
0(3) spin-ordered phase and spin-disordered phase. The spin-ordered phase is the Higgs phase of 
U(l) gauge dynamics, whereas the disordered phase is the confinement phase. We find a crossover 
in the confinement phase which separates dense and dilute regions of instantons. On the critical 
line, spin excitations are gapless, but the gauge-boson mass is nonvanishing. This indicates that a 
confinement phase is realized on the critical line. To confirm this point, we also study the noncompact 
version of the model. A possible realization of a deconfinement phase on the criticality is discussed 
for the CP JV +U(1) model with larger N. 



I. INTRODUCTION 

The concept of gauge theory plays an important role 
not only in elementary particle physics but also in 
condensed matter physics. In particular, in strongly- 
correlated electron systems and quantum spin systems 
at quantum phase transitions, unconventional low-energy 
excitations with "fractional" and/or "exotic" quantum 
numbers are expected to appear. In order to describe 
these excitations, gauge-theoretical approaches are some- 
times quite useful. 

For example, in the fractional quantum Hall effect 
(FQHE), each electron in a strong magnetic field sep- 
arates into a so-called composite fermion (CF) and even 
number of fictitious magnetic flux quanta. 1 Various ex- 
periments indicate that the CF's, instead of electrons, 
behave as quasiparticles. In a gauge theoretical approach 
to describing this phenomenon, a gauge field binding each 
CF and flux quanta is introduced 2 . When the dynamics 
of this gauge field is realized in a deconfinement phase, 
CF's and flux quanta are to dissociate each other, and 
the CF's become quasiparticle moving in an effective 
magnetic field made of the external magnetic field and 
the condensation of flux quanta. Then the integer QHE 
of CF's explains the FQHE. This mechanism of FQHE 
was called particle- flux separation (PFS), and the criti- 
cal temperature of confinement-deconfinement transition 
below which the CF picture holds was estimated 2 . 

In the high-T c superconductivity, a phenomenon called 
charge-spin separation (CSS) 3 was proposed to explain 
the anomalous behavior of the metallic phase of cuprates. 



In the CSS scenario, each electron in a cuprate is viewed 
as a composite of a holon and a spinon. A gauge- 
theoretical understanding of the CSS parallels the PFS 
as a deconfinement phenomenon of the gauge field bind- 
ing a holon and a spinon 4,5 . Then holons and spinons 
themselves behave as quasi-free particles. 

Recently, importance of quantum phase transitions 
(QPT) has been recognized, and several interesting pro- 
posals have been made 6 . One of the most tractable sys- 
tems of QPT is quantum spin models in low dimensions. 
By varying some parameters of a quantum spin Hamil- 
tonian, a phase transition from one phase with a certain 
order (e.g. the Neel order) to another phase with different 
order (e.g. the dimer) may occur. Then one is interested 
in what is the low-energy effective theory of the system 
that corresponds to the Ginzburg-Landau theory (GLT) 
of a second-order phase transition, and how it describes 
the phase transition. One may also ask the nature of the 
low-energy excitations appearing just on the criticality. 

For the conventional superconductivity, the GLT is 
well known and involves a collective scalar field describing 
Cooper pairs and a gauge field describing electromagnetic 
(EM) vector potential. The superconducting phase cor- 
responds to the Higgs phase of the EM gauge dynamics, 
in which gauge bosons (photons) acquire a mass via the 
Anderson-Higgs mechanism. The normal phase corre- 
ponds to the Coulomb phase with massless gauge bosons. 

For a certain class of quantum s = 1/2 spin models 
in two dimensions, it has been proposed that the corre- 
ponding "GLT" is a CP 1 model coupled with a dynamical 
U(l) gauge field 7 . This is quite natural because the path- 
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integral representation involves a set of coherent states of 
s = 1/2 quantum spin, which is conveniently labelled by 
a CP 1 (complex projective) variable, sometimes called 
Schwinger boson. The U(l) gauge field here appears as 
s = collective excitations to generate effective interac- 
tions between CP 1 spinons. 

Then it is an urgent matter to set up a concrete CP 1 
model coupled with a U(l) gauge field, and study its 
phase structure, low-energy excitations, etc. In this pa- 
per, we shall consider a three-dimensional (3D) lattice 
CP 1 model and study its properties by numerical meth- 
ods very intensively. The model contains two parameters 
ci and C2; ci controls the spin stiffness (the hopping pa- 
rameter of CP 1 bosons) and C2 is the inverse gauge cou- 
pling constant which controls fluctuations of the U(l) 
gauge field. This model may be viewed as an effective 
model of a quantum spin system on a 2D square lattice 
at zero temperature (T). It is also closely related to the 
t-J model of the cupratcs 4 ' 5 and the gauge-theory ap- 
proach to heavy fermion systems 8 . 

Let us briefly explain the relation between the present 
model and the t-J model of strongly correlated electrons. 
In the previous papers 9 , we studied the t-J model in the 
slave-fermion (SF) representation, in which spinons are 
described by CP 1 bosons and holons are described by 
fermions. At the half filling (no holons), the system be- 
comes the 2D s = 1/2 antiferromagnetic (AF) Heisen- 
berg spin model. At T = it is described by the 3D 
CP 1 model with infinite gauge coupling, which is just 
the present model at c-i = 0. The model exhibits a phase 
transition at C\ = C\ c between the 0(3) spin-ordered and 
disordered phases 4 . (See Sec.2B and Sec. 4 for more de- 
tails.) In lightly doped cases, holons can be integrated to 
"rcnormalize" the spin-stiffness constant c\ and destroy 
the long-range Neel order as the holon density increases 9 . 
For larger dopings, holons generates nontrivial interac- 
tions between the gauge bosons. Such a system may be 
mimicked to certain extent by the present model with a 
suitable amount of the C2-term [See Sec.2B(ii)], although 
one needs certainly more quantitative study on this point. 

The rest of the paper is organized as follows. In Sec. 2, 
we shall introduce the model and explain its relations 
to the classical 0(3) and 0(4) nonlinear sigma models, 
the GLT for quantum spin models, etc. In Sec. 3, we 
study the phase diagram by mean- field theory (MFT). 
In Sec. 4, we calculate the specific heat and present the 
global phase diagram of the model in the C2 — c\ plane. 
The model exhibits a second-order phase transition along 
the critical line c\ — C\ c (c2), which separates confinement 
phase and Higgs phase. In Sec. 5, we study spin-spin cor- 
relation functions. The results shows that the observed 
phase transition accompanies a spontaneous magnetiza- 
tion and the spin excitations are gapless on the critical 
line as well as in the spin-ordered phase. In Sec. 6, we 
calculate instanton (monopole) density in various places 
in the phase diagram. This supplies us with useful in- 
formation to understand the gauge dynamics. In Sec. 7, 
we measure the gauge-boson mass by calculating gauge- 



invarinat correlation functions of field strength. Two 
mechanisms of generating a gauge-boson mass are ex- 
plained; (i) condensation of instantons and (ii) Anderson- 
Higgs mechanism by condensation of the CP 1 variables. 
We also obtain the gauge-boson mass in X(f> 4 Higgs mod- 
els and compare the results with the previous calculations 
in order to check the present calculations. Section 8 is 
devoted to discussion. 

II. THE 3D CP 1 +U(1) LATTICE GAUGE MODEL 
A. Model 

Let us define the CP 1 +U(1) gauge model on the cubic 
lattice. We use x as the site index, and /j, = 0, 1, 2 as the 
direction index. We use \i also as the unit vector in the 
/i-th direction. On each site x we have a CP 1 variable 
z x , which is a two-component complex number, 

z x = ( ) , Z X 1,Z X 2 e C, (2.1) 

satisfying the CP 1 constraint, 

X>*»| 2 = 1 - ( 2 - 2 ) 
0=1 

On each link (x, x + 11) we have a U(l) gauge variable, 
U xfJ , = exp(i6 xfl ) [6 xfl e (-7T, +71-)]. 

The partition function Z of the model is written as 

Z = J [dU] J [dz] exp(-S), 

x.fj, x,/j, 

[dz] =Y[dz xl dz x2 S(\z xl \ 2 + \z x2 \ 2 - 1). (2.3) 

X 

The action S of the model is given by 

$ = 2" ^ ^ (^x+UfCtUxpZxa H~ H.C.^ 
x.fi.a 

~ ]T (u xv U x+ ^U x+ ^ v U Xft + -ll.c?j, (2.4) 

x,p,<v 

where c\ and c 2 are real parameters of the model. 

The action S is invariant under the following U(l) local 
(x-dependent) gauge transformation; 

Zxa > Z xa — GX^){iKx)z xa: 

U xli -> U' xli = cxp(iA x+l t)U xll exp(-iA x ). (2.5) 

A genuine CP 1 variable is characterized by (i) the con- 
straint (2.2) and (ii) the redundancy under the phase 
transformation z xa — ► cxp>{ik x )z xa . The gauge symme- 
try (2.5) assures us of this redundancy. 
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Besides the above local gauge symmetry, S has a sym- 
metry under the following global SU(2) transformation, 



Vz x ,VeSV{2). 



(2.6) 



This induces the global 0(3) spin rotation, which shall 
be explained in the following subsection B(ii). 



B. Special cases 

The present model reduces to some known models in 
certain limits of the parameters. 

(i) Pure gauge model (ci = 0) 

For ci — the model reduces to the 3D compact U(l) 
gauge model without matter couplings, where ci is the 
inverse of the gauge coupling constant <?, i.e., C2 oc g~ 2 . 
Polyakov 10 showed that this model stays always in the 
confinement phase regardless of the value of C2 due to 
the condensation of instantons (monopoles). We shall 
study the instanton effect and verify his result by numer- 
ical calculation in Sec. 6. 

(ii) 0(3) spin model (c2 = 0) 

For c 2 = the model reduces to a 0(3) nonlinear sigma 
model (classical Heisenberg spin model) with a modified 
action. Here one can perform the integration over C/ I(1 ex- 
actly because the plaquette coupling between U x ^s dis- 
appears. One obtains 



Z= f [d5]exp[]Tm/o(ciWi 
with the following measure, 

^ = nn*E sa xi-i). 



+ Sx+u ' S x 



)]> (2-7) 



(2.8) 



x i—1 i—1 



where Io(l) is the modified Bessel function. S x is a three- 
component 0(3) clasical spin vector, made of the CP 1 
variable, 



Zxcyzxi S x S x 1 , 



where a is the 2x2 Pauli spin matrices. We note that 
S x is a gauge-invariant object. A global 0(3) rotation of 
S x , which is a symmetry of the model, is induced by a 
global SU(2) transformation of z x , 



s x ^ns x ,neO(3), 

z x ^Vz x ,V eSU(2). 



(2.10) 



In Eq.(2.9) we have used the relation [dz]f({z x az x }) = 
[dS]f({S x }). Since ln/o(7) is a monotonically increasing 
function of 7, it is almost obvious that the model (2.7) be- 
longs to the same universality class as the standard 0(3) 



model with the action — J^2 S x +^ • S x . Then the model 
(2.7) should exhibit a second-order phase transition at a 
certain critical value c\ = ci c (0(3)), which separates the 
0(3) spin-ordered phase and the disordered phase. It is 
obvious from Eq.(2.7) that the parameter a controls the 
spin stiffness. 

For small C2, integration over the gauge field U xf i can 
be done perturbatively in powers of C2. The plaquette 
term in Eq.(2.4) generates four-spin couplings like (S x +ij,- 
S x ){S x +n+i/' S x + V ) in the leading order of C2, which prefer 
a ferromagnetic order for > 0. The higher-order terms 
of C2 become nonlocal and contain many spin variables. 

For any finite value of c 2 , integration over the gauge 
field U xfl in the partition function Z in (2.3) results in an 
effective spin model Ssm{S x ), 



[dS] exp(-S SM ) = / [dU] / [dz] exp(^S), (2.11) 



because "spin variable" S x is the only gauge-invariant 
object made of z xa . In later discussions, we call the 
off-diagonal long-range order of S x magnetization since 
it measures a spontaneous breakdown of the SU(2) or 
0(3) global symmetry (2.10). Later numerical calcula- 
tion shows that there exists a second-order phase transi- 
tion in the present model for which S x works as an order 
parameter. 

(iii) 0(4) spin model (c 2 =00) 
In the limit c 2 — > 00, fluctuations of U XI1 are totally sup- 
pressed and U x ,fi is restricted to pure gauge configura- 
tions, i.e., Uxfj, = cxp(— iAx+fj.) exp(iA x ). Then Z is re- 
expressed in terms of four-component 0(4) classical spin 
variables R xa {oi = 1, 2, 3, 4) G R, 

z x i — Rxl + iRx2, Z X 2 — Rx3 + iRxi, 
4 

Z X Z X — ^ ' Rxa = 1) 
a=l 

4 

[dz] = [dR] = n n dR ™ ■ *(E v (2.12) 



x a—1 



as 



( 2 - 9 ) Z = J [dR] cxp( Cl 



(2.13) 



This model is just the 3D nonlinear 0(4) sigma model 
which exibits also a second-order phase transition at 
ci = ci c (0(4)), which separates the 0(4) spin-ordered 
phase and the disordered phase. 



C. The nonuniform AF Heisenberg model 

Let us explain how the model Eq.(2.4) is related to the 
recently proposed effective field theories (GLT's) for the 
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quantum magnets. In Ref. 11 , we studied a phase transi- 
tion from the Neel state to the dimer state in the nonuni- 
form s = I AF Heisenberg model on a square lattice, 



HaF — ^ ' JxjSx • S x 



(2.14) 



where j is the spatial direction index (J = 1,2), S x is the 
quantum spin operator at site x and J X j is the nonuni- 
form exchange coupling. We rename the even lattice sites 
x = (o, i) where o denotes each odd site, and the index 
i = 1,2,3 and 4 specifies its four nearest-neighbor (NN) 
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FIG. 1. 2D square lattice; crosses are odd sites, and filled 
circles are even sites. Solid bonds indicate that their exchange 
couplings are stronger than those on dotted bonds. 



The exchange couplings J x j = J i may be position 
dependent, and in Ref. 11 we explicitly considered the fol- 
lowing dimer configuration, 



J oi 

AJ ni = 



J + AJoi, 

AJ ol = aJ, 
AJ ol = -aJ, 



i = 1 

i = 2,3,4 



(2.15) 



where < a < 1 is a parameter which interpolates be- 
tween the uniform Heisenberg model (a = 0) and the 
dimer model (a = 1). Its effective field theory is derived 
by the standard method 9 ' 12 in path-integral formalism 
by integrating over the odd-site spins. It takes a form of 
a CP 1 model in (2+l)-dimensions (the extra one dimen- 
sion is the direction along the imaginary time). Namely, 
the effective theory has the form of Eq.(2.4), where the 
effective parameter c\ is expressed in terms of a and the 
lattice spacing a as 



ci = 



1 I- a 2(2 + a) 



2y/2a 



1 



(2.16) 



The bare effective parameter C2 is vanishing in this effec- 
tive field theory. 

As explained in the case (ii) above, the model with 
C2 = exhibits a second-order phase transition at ci = 
ci c (0(3)). Through the correspondence (2.16), this tran- 
sition from the ordered state to the disordered state of 
the CP 1 model describes the transition of the nonuniform 
AF magnet from the Neel state to the dimer state 13 . In 
Ref. 11 we pointed out the possibility that the spinons 
z X a (a = 1) 2) appear as gapless excitations at the critical 



point. That is, a deconfinement phase like the Coulomb 
phase in the (3+l)-dimensions is realized there because 
of the appearance of long-range correlations among the 
gauge field. In this paper, among others, we shall reex- 
amine this point, i.e., investigate the possibility whether 
the Coulomb phase is realized on the critical line of the 
model (2.4). 



D. The uniform AF Heisenberg model with ring 
exchange 

In a framework of the uniform AF Heisenberg model 
with multiple spin ring exchange, Senthil et al. 7 argued 
that a phase transition from the Neel state to valence 
bond solid (VBS) takes place. They suppose that the 
system in the Neel state but near the critical point is 
described by a CP 1 field coupled to a noncompact U(l) 
gauge field, where nonperturbative contributions from in- 
stantons are totally suppressed. An effective field theory 
in the continuum space-time in this region was given by 7 



"deconfine 



\(d^ ~ iA^)S a \ 2 + s\2 a \ 2 + u(\S a \ 2 ) 
9 



(2.17) 



where z a (a = 1,2) are unconstrained complex scalar 
fields coming from the CP 1 spin variables, is the gauge 
potential, and s and u are parameters. For s < s c , z a 
are condensed (z a ) ^ which corresponds to the long- 
range Neel order. The gauge dynamics is in the Higgs 
phase and the gauge bosons acquire a mass gap \ (z)\. At 
the criticality s = s c , (z a ) = 0, whereas the effect of 
instantons is negligible because of the cancellation mech- 
anism via the Berry phase term 14 . On the other hand, in 
the VBS state, the instatons are relevant and proliferate 
to generate the confinement phase. The proposed phase 
diagram 7 of the AF Heisenberg model is shown in Fig. 2. 
We shall compare it with that of the present model in 
Sec.7C. 



VBS 



Neel 




U(l) spin 
liquid 





FIG. 2. Renormalization-group flows of the effective gauge 
theory of the AF Heisenberg model with the ring exchange 7 . 
The horizontal axis refers to the parameter s, and the vertical 
line to the instanton fugacity /. 
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III. MEAN FIELD THEORY 



Before going into detailed numerical studies of the 
model, let us apply a MFT to the model to obtain a 
rough phase diagram. 

The MFT below is based upon the variational 
principle 15 . We start with the variational action So and 
use the following relations; 

Z = J [dz][dU] exp(-Sb) = exp(-Fo), 
F=-lnZ <F V = F + (S-So)o, 
(O) = V J [dz][dU] O cxp(-S* ). (3.1) 

From this Jensen-Peierls inequality, we adjust the varia- 
tional parameters contained in So optimally so that F v 
is minimized. 

For the trial action So, we assume the translational 
invariance and consider the following sum of single-site 
and single-link energies; 

x a 

where W and H are complex variational parameters. 
Then we obtain the following free energy per site, /„ = 
F v /N, where N is the total number of the 3D lattice sites 
(we present the formulae for d-dimensional lattice); 

h=\H\, u = \W\, 

fv = -dln(I (u)) - \n(I (h) - I 2 {h)) - Cl dm 2 p 
diop + hm, 



d(d-l) 4 
C 2 n P 



,^ Hz v ,2,1/2 = j h(h)-h(h) 
l ^ IWo1 ' ~ 2I (h)-I 2 (h)' 



P = \(U xli )o\ 



AM 



(3.3) 



where I n (7) (n is an integer) is the modified Bessel func- 
tion, 

r2iK 

!n{l) = I ^- exp(7 cos 9 + ine). (3.4) 
Jo 27T 

The stationary conditions for f v w.r.t. u>, h read 

us = cim 2 + 2c 2 (d - l)p 3 + 6c 3 (d - l)m 2 p 2 , 

h = 2dc\mp + 4c 3 d(d - l)mp 3 . (3.5) 

The MFT equations (3.3)-(3.5) for d = 3 predict the 
existence of the three phases characterized as follows; 



phase 


P = (U Xf i) 


m = (z x ) 


Higgs 




^0 


Coulomb 







Confinement 









(3.6) 



The naming of each phase should be clear. We note that 
the fourth combination p = and to 7^ is missing. 



Higgs 



Confinement 



Coulomb 



FIG. 3. Phase diagram of the 3D CP 1 +U(1) model in 
MFT. There are three (Higgs, confinement, Coulomb) phases. 



In Fig. 3, we plot the phase boundaries obtained from 
Eqs.(3.3)-(3.5). The phase boundary of the MFT be- 
tween the Higgs and Coulomb phases is of second order, 
whereas other two boundaries, the Higgs-confinement 
and confinement-Coulomb, are of first order. Across the 
second-order transition, p and to vary continuously, while 
across the first-order transitions, p and/or to change dis- 
continuously with finite jumps of Ap and/or Am. For 
the Higgs-confinement transition, Ap ^ and Am 7^ 0, 
and for the confinement-Coulomb transition, Ap ^ and 
Am = since to = in both phases. 

The present MFT for the CP 1 + 17(1) model may have 
some inappropriate points: 

(i) As explained above, for the pure gauge case c\ = 0, 
the confinement phase should survive up to c 2 — > 00 and 
the Coulomb phase in the MFT should disappear. In fact 
we verify the nonexistence of the Coulomb phase in the 
following section. However we shall see some crossover 
phenomenon close to this ficticious transition line. See 
later discussion. 

(ii) At C2 = 0, the system reduces to the 0(3) spin 
model (2.7) having a second-order transition at c\ = 
c lc (0(3)). Thus the first-order transition of the MFT at 
c-2 = is incorrect. This MFT result should be modified 
by Monte Carlo simulations in the following section. 

On the other hand, in the limit c 2 — > 00, the MFT 
correctly predicts the order of transition, i.e., the second 
order one of the 0(4) spin model. 



IV. PHASE STRUCTURE 

To study the phase structure of the model numerically, 
we performed Monte-Carlo simulations with Metropolis 
algorithm for the 3D cubic lattice of the size up to N — 
L 3 = 30 3 with the periodic boundary condition (PBC). 

We first calculated thermodynamic quantities per site 
like the internal energy U and the specific heat C, 



U=(S)/N, 
C=((S-(S)) 2 )/N 7 



(4.1) 



5 



as functions of the parameters c\ and c 2 . The specific 
heat C has a peak as a function of c\ for each fixed 
c 2 . See Fig. 4 for c 2 = and Fig. 5 for c 2 — 2.0. These 
peaks develop as the lattice size increases as N = 6 3 , 10 3 , 
20 3 and 30 3 for the c 2 = case and N = 8 3 , 10 3 and 
20 3 for the c 2 = 2.0, respectively. This behavior indi- 
cates second-order phase transitions. Fig. 4 shows that 
the 0(3) spin model of Eq.(2.7) has a phase transition at 
c lc (0(3)) = 2.85.... 



CI 

2 2.25 2.5 2.75 3 3.25 3.5 
FIG. 4. System-size dependence of the specific heat for 
c 2 = 0. 



C N (e) = N^^N^e), 



(4.3) 
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where Cjy(e) is the specific heat in the system of size 
TV and 4>{x) is a scaling function. We show the scaling 
functions <f> m Fig. 6 which are determined from the data 
in Figs. 4 and 5. The parameters of the FSS are = 
2.87, v = 1.50, a = 0.20 for c 2 = and c loo = 1.09, v = 
1.17, a — 0.17 for c 2 = 2.0, respectively. We think that 
the finite-size scaling hypothesis (FSSH) holds quite well 
for the c 2 = case. For the c 2 = 2.0 case, data exhibit 
the FSS in smaller parameter region close to the peaks 
of C than that of the c 2 = case. We think that it is 
due to the existence of another length scale besides the 
spin correlation length £, i.e., the correlation length of 
the gauge boson. See the discussion in Sec. 7. 

Next, in Fig. 7, we plot C at ci = as a function of 
c 2 . It has certainly a peak at c 2 = 1.4..., which, how- 
ever, does not develop as the lattice size increases. Thus 
this peak does not indicates any phase transition along 
ci = 0. This nondeveloping peak at (ci = 0, c 2 = 1.4...) 
continues to the region c\ > 0. This line of the peaks may 
exhibit some crossover. The physical meaning of this line 
will be discussed after calculating instanton density in 
Sec.6. 
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FIG. 5. System-size dependence of the specific heat for 
c 2 = 2.0. 



It is interesting to see if the data of the specific heat 
exhibit the finite-size scaling (FSS) behavior. Let us in- 
troduce a parameter e as 



where cioo is the critical coupling at the infinite volume 
limit. The correlation function £ behaves as £ cx e~ u 
with a critical exponent v. It is also expected that the 
peak of C diverges as oc t~ a as N — > oo with an- 
other exponent a. Then for sufficiently large systems, we 
expect 



FIG. 6. Scaling function <j>(x) obtained from the data of 
the C2 = (left) and C2 = 2.0 (right) cases. 
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12 3 4 

FIG. 7. System-size dependence of the specific heat peak 
for ci = 0. No developments of peaks are observed. 



In Fig. 8 we present the phase diagram in the ci — c\ 
plane obtained via these specific-heat measurements. We 
also show the locations of the nondeveloping peaks as 
crossovers. 



♦ :Second order phase transition 
x :Crossover 



Higgs 




confinement 



C2 



1 2 3 4 5 

FIG. 8. The phase structure of CP 1 + (7(1) model in the 
C1-C2 plane obtained by the measurement of the specific heat. 



Cs(t) exhibits a nonvanishing off-diagonal long-range or- 
der (magnetization) M, M = lim^oo Cs {t) , where 
t — ► 00 implies t = L/2 for a finite lattice of size N = L 3 . 
Cs(t) approaches to M 2 algebracally, 



c s (t) = M 2 + Ar { - 1+r ''\ 

rf = 0.0 - 0.1 for L = 20. 



(5.5) 



In contrast, at a lower c\ = 2.5 [< ci c (0(3))], M = and 
Cs(t) decays exponentially, 



C s (t) = A'expi-j't). 



(5.6) 



Thus the phase transition observed in the previous sec- 
tion corresponds to a spontaneous symmetry breaking 
of the SU(2) spin symmetry (2.10) of z x , say (z x ) — 
(M 1 / 2 , 0)*. To confirm this observation, we plot the mag- 
netization M for C2 = in Fig. 10. It exhibits a typical 
behavior of second-order phase transition. 
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FIG. 9. Spin-spin correlations for C2 = 0. For c\ — 2.5, an 
exponential decay A' exp -7 ' fits it better. For ci = 2.8 and 
3.4, an algebraic decay M 2 + At~ a fits it better. 



V. SPIN EXCITATIONS 

In this section we study the spin-spin correlation func- 
tions and the energy gap of spin excitations. The corre- 
lator of 0(3) spins, S x of Eq.(2.9), is defined as 



M 



Cs(t) = {S x ■ S x +tij)- 



(5.4) 



Cs(t) is a function of the distance t but does not depend 
on the reference point x nor the direction fi due to the 
PBC on the cubic lattice. 

Let us first consider the case C2 = 0. As explained be- 
fore, this case corresponds to the nearest-neighbor spin 
model (2.7). In Fig. 9 we plot Cs(i) for the lattice size 
20 3 . Cs(t) changes its behavior around c\ = ci c (0(3)) = 
2.8..., the critical point determined by the location of spe- 
cific heat of Fig.4. At a higher a = 3.4 [> ci c (0(3))], 
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FIG. 10. Spontaneous magnetization for C2 = 0. The result 
shows a typical behavior of second-order phase transition. 



In Fig. 11 we also show the magnetization M at various 
values of C2- M starts to develop at the critical points 
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observed by the specific heat measurement. This means 
that M is an order parameter for the phase transition in 
the present model. 
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FIG. 11. Spontaneous magnetization for various C2. The 
result shows a typical behavior of second-order phase transi- 
tion. 



As we mentioned in Sec. 2, the model reduces to the 
0(4) spin model in the limit c 2 — > oo. We studied the 
phase structure of this 0(4) model, and found that there 
exists a second-order order-disorder phase transition as 
in the pure CP 1 model. The critical coupling is estimated 
as ci c (0(4)) = 1.0... which is fairly close to the critical 
value of c\ for sufficiently large c 2 . We also investigate 
the spin correlation function in the ordered phase. The 
data can be fitted well in the form (see Fig. 12), 



C s (t) = M 2 + br ( - 1+r >"\ 
rf = 0.1 - 0.3 for L = 20. 



(5.7) 



This result indicates that smooth transfer in the behavior 
of Cs(t) occurs from the 0(3) spin model at c 2 ~ to 
the 0(4) spin model at c 2 — oo. 



C s (t) 




0.175 



t 



1 2.5 5 7.5 10 12.5 15 17.5 
FIG. 12. Spin correlation function in the 0(4) model at 
ci = 1.60 and C2 = oo. The system size L = 20, and the 
exponent is estimated as rf = 0.26. 



Next let us study the mass (energy gap) Ms of spin 
excitations. To measure Ms with good precision, we in- 



troduce a Fourier transform of the spin variables S x in 
the 2D plane (xi,X2) 



(5.8) 



In the continuum limit, one can readily find that the 
correlator of S(xa ; pi , P2) behaves as 

(S(x 3 ;pi,p 2 ) ■ S(0;-pi,-p 2 )) = / dp, 



■p>+Ml 



oc e -VPl+P 2 2+ M s^ : (5.9) 

where p 2 — Yli=i2 3Pi an d Ms is the energy gap of 
the lowest spin excitations. Below we fit the measured 
correlation function on the lattice by using this form to 
calculate Ms- Let us define 

D s(t) = ^2(S(x 3 + t;p x ,p2) ■ S(x 3 ; -pi, -p 2 ))- 



(5.10) 



The PBC restricts the momentum as pi = 2irrii/L, m = 
0, 1,2, — 1. In the simulation we choose nonvanish- 
ing smallest values, p\ — P2 — 2tt/L in order to fit the re- 
sult in an exponentially-decaying form even for Ms = 0. 
We note that, from the above definition of Ms, it is pos- 
sible that Mg takes a negative value as long as the ex- 
ponent p\ +P 2 , + Mj < 2(2n/L) 2 . We shall comment on 
this point later on. 

In Fig. 13, we plot M s for TV = 20 3 . For a fixed c 2 , M s 
is a smooth function of c\ and vanishes in the spin ordered 
phase c\ > c\ c as predicted by the Nambu-Goldstone 
theorem. We checked the result of Fig. 13 by repeating 
calculation with the choice {p\,P2) — (2ir/L i 0). We shall 
use the same techniques for calculating the gauge-boson 
mass in the following section. 
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FIG. 13. Energy gap Ms of 0(3) spin excitations for 
C2 = 0.0, 1.2 and 2.0. The errors become larger for smaller 
ci because of more rapid fall off of Ds (t) . 
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VI. INSTANTONS 
A. Definition of instanton density 

In this section, we calculate the instanton density Q x at 
x in various set of (c2,c\). Measurement of Q x provides 
us with important information about fluctuations of the 
gauge field, though it is rather difficult to determine a 
location of phase transition by itself. 

We follow the definition of integer instanton charge 
by DeGrand and Toussaint 16 . First, let us consider the 
magnetic flux Q Xt)lv penetrating plaquette (x, x + fi, x + 
(X + V, x + v), 

(-4tt < e XjM „ < 4tt). (6.1) 

We decompose Q x ,nv into its integer part 2nn XtfiI/ 
{n x ,fiv is an integer) and the remaining part Q XtfiU = 
®x,nv (mod 2ir) uniquely, 



flu 



e 



X,fJ,V 1 



(-7T < Q Xtflv < 7T). (6.2) 



Physically, n x ,nv describes the Dirac string whereas Q x ^ v 
describes the fluctuations around it. The quantized in- 
stanton charge Q x at the cube around the site x = 
- 1 - of the dual lattice is defined as 



2 1 2 1 2 



Qx ~ ^ ^■fJ,up{j^X-\-jl,up TlxjUp) 



^ ^ ^ £ p,v p{® x+ p,v p ®x,vp 



(6.3) 



in Fig. 14]. The result indicates that the instanton den- 
sity decreases very rapidly from C2 ~ 1.5 as ci increases. 
Therefore, we can identify the crossover found by the 
specific heat C in Sec. 4 as a crossover from the region of 
dense instantons to that of dilute instantons. In fact this 
was first observed in Ref. 17 . 
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FIG. 14. Instanton distribution is measured at various 
places (a) ~ (1) in the phase diagram in the C2 — ci plane. 



where e^up is the complete antisymmetric tensor. 
Q x measures the total flux emanating from the 
monopole (instanton) sitting at x. Roughly speaking, Q x 
measures the strength of nonperturbative gauge config- 
urations. Polyakov showed that a condensation of these 
instantons drives the system at c\ = into a confine- 
ment phase, which is characterized by strong and large 
fluctuations of gauge field 8 xfl at long distances. 



B. Results 



In this subsection, we shall show "snapshots" of Q x , 
typical configurations of Q x taken in the process of MC 
updates after sufficient thcrmalization. In Fig. 14, we first 
show the locations of pairs of (02, ci) in the phase diagram 
at which snapshots are taken. In Fig. 15, we present the 
snapshots at these locations for the lattice N — 16 3 . In 
Fig. 16, we show Q{n), the density of instantons with the 
charge n(> 0), Q(n) — J2 X &n,\Q x \/N. The length of each 
hexahedron shows the magnitude of instanton charges 
there. 

Let us first see the behavior of instantons in the pure 
3D U(l) gauge system with c\ = [(a),(b),(c) and (d) 
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(i):cl=0.5, c2=2.0 (j):cl=2.0, c2=2.0 (k)xl=1.0, c2=1.0 (l):cl=2.5, c2=1.0 

Q(1)=Q(-1)=10 Q(1)=Q(-1)=0 Q(1)=Q(-1)=374 Q(1)=Q(-1)=20 

FIG. 15. Snapshots of spatial distribution Q x of instantons. Dark boxes indicate the locations of Q x > whereas light boxes 
indicate Q x < 0. Q(n) = Y] x 5 n ,|Q x | is the density of instantons with the charge n. 
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Q(n)/N Q(n)/N 




FIG. 16. Distribution of the instanton charge density Q(n) 
at the locations (e) — (k) in the phase diagram. Q(n > 3) = 0. 



More systematically, we calculate the average density, 
p = (J2 X \Qx\}/N for ci = as a function c%. See Fig. 17. 
We fit the result with exp(— Ic-i) as expected by the 
dilute-gas approximation 10 , where £ (~ 4.35), whereas 
its theoretical value is calculated as 5.06 10 . The fitting is 
quite satisfactory for C2 > 1.5 as expected. 



p 




FIG. 17. Average of instanton density p = (Y] \Qx\)/N 
for ci = as a function of C2- The solid curve is a fit in the 
form of dilute-gas approximation oc exp(— £02) with £ ~ 4.35. 

Let us turn to the pure CP 1 model at C2 = 0. In Fig. 15, 
the snapshots ((e) ~ (h)) represent this case. The ci-term 



of the action shows that a spin order, i.e., a condensa- 
tion of z x , suppresses the fluctuations of U xfJi . Due to 
this correlation z x+tl z x «-> U xfi , instanton configurations 
in (2+l)-dimcnsions correspond to creations and/or de- 
structions of so-called skyrmions configurations of 0(3) 
spins in the 2D plane. Therefore instantons are sup- 
pressed in the spin-ordered state. The results in Fig. 15 
are consistent with the above observation as the instan- 
ton density decreases as c\ increases. However, we note 
that even in the spin-ordered phase at c\ = 3.5 there are 
a number of surviving instantons, most of which make 
instanton-anti-instanton dipoles in the nearest-neighbor 
pairs. 

Finally, let us take a look at the snapshots obtained in 
the region close to the critical line; the snapshots at (i) ~ 
(1) in Fig. 15. We see that the instanton density increases 
as ci and/or C2 decrease. Even in the spin-ordered phase, 
the instanton density is finite for small C2. We do not 
think that this result contradicts a natural belief that 
the Higgs phase is realized in the spin-ordered phase. 
For the confinement phase to appear, finite instanton- 
density is not sufficient, but dissociations of instanton- 
anti-instanton dipoles are necessary. 

The observation of instantons in the present section is 
quite helpful for understanding the behavior of gauge- 
boson mass, which is calculated in the following section. 

VII. MASS OF GAUGE BOSONS 

In this section we study the gauge-invariant mass Mq 
of gauge bosons. The present model contains two inde- 
pendent dynamical variables, i.e., the CP 1 field and the 
gauge field, and therefore there exist two independent 
length scales. In Sec. 5, we studied the correlation length 
of the CP 1 spin sector. In this section, we investigate 
that of the gauge sector. 

We explain the definition of Mq in Sec. 7.1, and show 
the numerical result of Mq in Sec. 7. 2. In Sec. 7.3 we 
address the problem of whether the system is in the con- 
finement phase or in Coulomb phase just on the critical 
points c\ c . In Sec. 7. 4, we study the case of noncompact 
gauge field to support our argument for the problem in 
Sec.7.3. 



A. Definition of gauge-boson mass Mq 

In this section, we study the gauge-invariant mass of 
gauge bosons, Mq, by using the same techniques that we 
used to measure the spin gap Ms in the previous sec- 
tion. To define Mq we first introduce a gauge-invariant 
operator 0(x), 
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= 51 e^v sin(-0 x „ - O x+Vtll + 8 x +n,v + Ox?)- (7-1) s Ka = - 

Then we intorduce a Fourier transform as before, 

6(13) = 0{x)e ipiXl+lp2X2 . (7.2) 



to define the gauge correlation function as 

D G® = ^3 £ (0(0:3)0(13+0 



(7.3) 



In the continuum, Dc(t) is expected to behave as 



De{t) 



dp-. 



p 2 



(7.4) 



We determine Mq by fitting the data in this expo- 
nential form (7.4). For practical calculations, we set 
Pi — P2 — 2?t/L as before. Here we note that we are 
assuming the ordinary form of the gauge-boson propa- 
gator for a gauge system coupled with massive matter 
fields. On the critical line, the gauge-boson propagator 
may be modified due to the appearance of massless field 
z xa . See later discussion on this point. 

In order to verify that the above methods give correct 
results, we first apply them to the 3D U(l) lattice gauge- 
Higgs models for which the gauge-boson mass has been 
already calculated by using other methods. In Ref. 18 the 
action S of the model is parameterized as follows; 



** = -|£IItf-f£(* 



K 



■£[I*»I 2 + A (I* 



X.fi 

2 



x+nU xfl $ x +H.c.) 
I) 2 ] , (7-5) 



where & x is the complex Higgs field. By changing the 
value of K, the system undergoes a Higgs-confinement 
phase transition. A gauge-boson mass was calculated 
by studying the genuine gauge-boson propagator with a 
specific gauge-fixing condition. We calculated the gauge 
boson mass for the same parameter region as in Ref. 18 
but by using Da{t) of Eq.(7.3). In this region, the gauge 
coupling constant has a moderate value (3 = 2 and A = 
0.020. In Fig. 18 we show our results of the gauge-boson 
mass, which are in good agreement with those of Ref. 18 . 
We should notice that the gauge boson mass is always 
nonvanishing even at the critical point, though one might 
expect the behavior Mq ~ \K — K c \ where K c is the 
value at the critical point. We shall comment on this 
point after showing the results of the present CP 1 +U(1) 
model. 

On the other hand, in Ref. 19 , it was reported that the 
gauge-boson mass vanishes in a certain parameter region. 
In Ref. 19 , the Higgs part of the action is parameterized 
in a complicated way as 



§ En u - f e (**+^*« + h - c -) 



K 
~2 



Ei^i 2 [ e 



X,fX 

y 3.17(1 + 2a;) 
K 2 2^K 



(-4 + 8x - 8ir 2 )(log6/3 + 0.09) + 25.5 + 4.6a; 



16tt 2 /? 2 



xK 2 



£i*« 



(7.6) 



We also calculated the gauge-boson mass of Eq.(7.6). In 
Fig. 19 we present our result of Mq for x = 2, (3 = 4 and 
K = 5, the same parameters used in Ref. 19 . The result 
agrees with that of Ref. 19 . For positive y, Mq reduces to 
around zero. This phenomenon may be understood natu- 
rally as follows; For large y the mass term of $ a becomes 
large and <i> x fluctuates weakly around zero. Then one 
can treat the hopping K term as a small perturbation. 
The main term is the first /3-term whose coefficient is 
large. Thus the fluctuations of gauge bosons are strongly 
suppressed, so one can expand as U x/Jl ~ l+iA xtl . The re- 
sulting action cx & x ^ lL/ describes just free massless gauge 
bosons. See also the discussion in the following subsec- 
tion. 
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FIG. 18. Gauge-boson mass of Eq.(7.5) for (3 = 2 and 
A = 0.020, the same parameters as in Ref. 18 . 
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FIG. 19. Gauge-boson mass of Eq.(7.6) for x = 2,13 = 4 
and K = 5, the same parameters as in Ref. 19 . Here and 
in the forthc oming figures of Mg, Mq denotes the product 
sgn(Mc) • y/\M%\ with M% defined by Eq.(7.4) in order to 
distinguish the case of negative Mg clearly. 



B. The case of compact U(l) gauge field 

We measured the gauge-boson mass for the 16 3 lattice. 
So the distance t in D G (t) of (7.4) takes t — l,---,8 
due to the PBC. The value of M G for each point of the 
C2 — c\ plane is calculaed as the average of 10 samples, 
each sample of which is obtained after 2 x 10 5 ~4x 10 5 
sweeps. The error of M G is evaluated by the standard 
deviation of these 10 samples. 

There are two mechanisms to generate nonvanishing 
Mq'- (i) Condensation of instantons and (ii) Anderson- 
Higgs mechanism. We shall see them first by studying 
Mq{c2) as a function of c 2 for a fixed Cx, and then by 
studying M G (c\) for a fixed 

In Fig. 20 we present M G (c<i) for c\ = 0,0.8 and 1.5. 
Let us first see the pure gauge case c\ = in some de- 
tail. As C2 increases from zero, Mq decreases and be- 
comes almost zero at c 2 ~ 2.5, the value at which the 
instanton density almost vanishes. Thus we find a pos- 
itive correlation between the gauge-boson mass and the 
instanton density. Actually, in the dilute-gas approxima- 
tion, Mq is estimated as Mq <x p cx exp(— Ic-i) 1 ® . In 
fact, for p ~ 0, one may forget the nonperturbative (in- 
stanton) effect and expands as U x ^ ~ 1 + i0 x ^ in the 
action. The action of small fluctuations 9 xfJj becomes 
S ~ Ylxuv®x)iu' wm ch describe a free massless gauge 
boson as explained before in the case of Higgs model 19 . 
On the other hand, for the system of dense instantons, 
the gauge field Q x>1 has short-range, i.e., massive correla- 
tions because of their wild fluctuations. In short, as the 
C2 term controls the density of instantons, i.e., wild fluc- 
tuations of 9 Xfl , Mq{c2) decreases as c 2 increases. Fig. 20 
shows that this tendency survives for c\ = 0.8 and 1.5. 
This is expected from Fig. 15 of instanton density. We 
remark here that the vanishing gauge-boson mass in the 
numerical calculations in finite systems does not guaran- 
tee that the system is in the deconfinement phase. The 
measurement of the specific heat indicates that the sys- 
tem is in the confinement phase even for large c 2 if ci 
is below ci c (c 2 ). We think that the above calculation 
Mq ~ in the confinement phase is a finite-size effect 
and if we calculate Mg in very large systems, we obtain 
Mg > in the confinement phase. 

In order to verify the above expectation, let us take 
a look at Fig. 20 in more detail. One sees that M G — 
-(0.1) 2 - -(0.15) 2 for c 2 > 2.5 with c x = and ci = 0.8. 
As mentioned below Eq.(5.10), these negative values of 
M G are possible in the present definition of M G . We 
think that this unphysical result stems from the finite 
size of the system as mentioned above. Similar results are 



reported in Ref. 20 for the D = 4 compact scalar QED. 
In order to see the finite-size effect on the gauge boson 
mass, we plot data of D G (t) m t — hiDG(t) plane in 
Fig. 21. From this plot, it is obvious that the PBC affects 
the behavior of D G (t) for t ~ L/2 where L is the system 
size. Then one may discard the data of t ~ L/2 and re- 
estimate the gauge-boson mass to obtain M G = — (0.08) 2 . 
To get a definite conclusion for the negative gauge boson 
mass, calculation in larger systems is required. In this 
direction, we calculated M G in a larger system 24 3 and 
obtained a reduced value M G ~ — (0.03) 2 . Therefore we 
expect that we obtain M G > by calculation in very 
large systems, and in particular M G > in the confine- 
ment phase as it is generally expected. 
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FIG. 20. Gauge-boson mass Ma as a function of C2 for 
ci = 0.0, 0.8, 1.5. The data of c\ = 0.0 and 0.8 almost overlap. 
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FIG. 21. Correlation function Dg{€) of gauge boson vs t at 
(ci = 0.0, C2 = 3.0) for L — 16. (a) the exponential curve fit- 
ting the data with t = dropped, (b) Log plot lnL>G(t). The 
solid straight line fits all the data, while the dashed straight 
line fits the data with t — 6, 7, 8 dropped. 

On the other hand, the data for c\ — 1.5 is interesting 
because the horizontal line with c± = 1.5 in the c 2 — c\ 
plane intersects the critical line as c 2 is increased. The in- 
tersection point is estimated as [c\ = 1.5, c 2 = 1.2) from 
the phase diagram Fig. 8. As seen in Fig. 20, the gauge 
boson mass decreases very rapidly until ci ~ 1.3 and be- 
comes almost constant for C2 > 1.4. This rapid decrease 
of Mq in the spin-disordered phase comes from the de- 
crease of the instanton density as we measured in the pre- 
vious section. The finite gauge-boson mass in the spin- 
ordered phase can be interpreted as due to the Anderson- 
Higgs mechanism because Mq stays almost constant for 
varying c^. We shall explain this mechanism in more 
detail for Mq{ci) below. 

Next, let us study Mq{c,\) for a fixed c%. In Figs. 22-26, 
we present Mq(ci). For all values of C2, as ci increases 
from zero, Mq{c,\) first decreases in the spin-disordered 
(confinement) phase, < ci < C\ c {c2), to reach the mini- 
mum at the critical point c\ — C\ c {c2). Then it increases 
in the spin-ordered (Higgs) phase ci c (c 2 ) < c\. The first 
decrease is due to the suppression of instantons at larger 



ci as we observed in the previous section. The situation 
parallels the pure gauge system ci = if the axes c\ and 
C2 are interchanged. 

The successive increase in the ordered phase is due to 
the Anderson-Higgs mechanism. In the ordered phase, 
(z x ) starts to develop at c\ = c± c as Fig. 10 shows. Then 
the c\ term of the action supplies the following mass 
term; 

Su = \M%Y^ e l^ M% = Cl {z x )\ (7.7) 

where we expanded U xfi ~ 1 + i0 xfi by assuming small 
gauge-field fluctuations 8 xfl . This expression Mq at a 
fixed C2 certainly increases as ci (> c\ c ) increases. 
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FIG. 23. Magnified figure of Fig. 22 near the critical line. 

An important point is that the gauge boson mass is 
nonvanishing on the critical line except for relatively 
large values of c 2 . To confirm this point we present the 
size-dependence of Mq in Figs. 24-26. They show that 
Mq slightly increases for a larger lattice, which denies 
the possibility that Min[Mc(ci)] — > as L — > oo. This 
point is in contradiction to the expression (7.7), where 
(z x ) 2 = at ci = cic implies Mg(ci c ) = 0. The nonper- 
turbative effects of U x ^ are crucial to explain Mq ^ 0. 
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In fact, we observe Mq at the criticality increases as C2 
decreases. This behavior is consistent with the result of 
the instanton density at c\ = c\ c , which is nonvanishing 
and increases as C2 decreases. This suggests that on the 
critical line, the confinement phase in which instantons 
survive, rather than the Higgs or Coulomb phases, is re- 
alized. In later subsection, we shall confirm this point by 
studying the noncompact version of the present model. 
Here we note that the fact Mg(ci c ) ^ does not con- 
tradict the general property of a second-order transition 
that the correlation length diverges at a critical point. In 
fact, we have observed in Fig. 13 that the spin gap Ms 
vanishes at c = c\ c . Thus the spin correlation length di- 
verges on the critical line, which is sufficient to generate 
singularities in thermodynamic quantities, although the 
gauge correlation length is finite there. 
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FIG. 26. Size dependence of the gauge-boson mass with 
c 2 = 3.0. 
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FIG. 25. Size dependence of the gauge-boson mass with 
c 2 = 2.0. 



C. Confinement phase vs Coulomb phase on the 
critical line 

Physical properties of a quantum system just on its 
critical point are interesting. For example in the previous 
paper 11 , by a gauge-theoretical study of a quantum spin 
system having a Neel-dimer (Higgs-confinement) transi- 
tion, we argued that the Coulomb- like phase is realized on 
the critical point. This is because there appears massless 
relativistic spinon z xa on the critical point and it gen- 
erates nonlocal correlations for the gauge field. Typical 
form of the effective gauge model is given as 



(7. 



where C denotes an arbitrary closed loop on the lattice 
and \C\ is its length. g e is the effective gauge coupling 
constant and 7 is a parameter of the hopping expansion. 
At the critical value 7 = 7c, which corresponds to rela- 
tivistic massless matter fields, S e s diverges for the uni- 
form configuration U x ,p — 1. Recently we studied models 
closely related to Eq.(7.8) and showed that there exists a 
dcconfincmcnt phase transition at 7 = 7 C and sufficiently 
large g e 21 ■ In the deconfinement phase, the renormalized 
gauge boson propagator behaves at long distance p ~ 
as follows; 



{A,{-p)A v {p)) 



1 



(7.9) 



which gives potential energy between charges at distance 
r as 1/r. Then we should evaluate the following Fourier 
integral instead of Eq.(7.4), 



F(p\+pl-t) = / al P3 ^ 
J VP 



D 2 +M 2 



(7.10) 



We numerical evaluated the above integral (7.10) and 
found 



F(e;t)~e- a,r/e+MH , a Q ~ 1.287. 



(7.11) 
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One may expect that the above change of the fitting 
function F(£; t) may explain the finiteness of the gauge- 
boson mass at the criticality. However this is not the 
case. Let us assume M G = on the critical line in the 
present model. Then the correlation function of 0(2:3) 
obtained by the propagator (7.9), D' G (t), behaves as 



D' G (t)\ MG=0 ~ e- Q °VpFpi 



(7.12) 



From Eq.(7.12) it is obvious that the function 
D' G (t)\Mo=0 is invariant under a scale transformation; 
(t,p) — > (£t,p/£) where £ is a patameter of the scale 
transformation. This means that, if Mq = and the 
"Coulomb phase" is realized at the critical line, the 
damping factor of D G {t) stays constant along the critical 
line. However in Fig. 22, the value of M G at minimum 
(which gives the damping factor) clearly increases as c 2 
decreases. This implies that there exists a physical scale 
of the mass dimension (i.e. the gauge-boson mass Mq) in 
the gauge-boson sector at the criticlity, which supports 
the conclusion in the previous section that Mq ^ on 
the critical line. 

In Fig. 27 we present the suggested flow diagram of the 
renormalization group (RG) in the present model. The 
RG flow on the critical line moves toward the unstable 
fixed point A in the phase diagram. The two flow di- 
agrams Fig. 27 and Fig. 2 are quite similar; The unsta- 
ble fixed point (s = s c , / = 0) in Fig. 2 corresponds 
to the point A in Fig. 27, and similarly the Neel point 
to the point (02 = 0, c\ = 00), the VBS to the point 
B(c2 = ci = 0), and the U(l) spin liquid to the point 
(02 = 00, ci = 0). More comments on this point will be 
given in Sec. 8. 




FIG. 27. Renormalization-group flow in the phase diagram. 
The thick line is the critical line. A is an unstable fixed point, 
and B is a stable fixed point where density of instantons di- 
verges. 



D. The case of noncompact U(l) gauge field 

For the U(l) gauge-Higgs model on a 3D lattice, ex- 
actly the same behavior of the gauge-boson mass as ours 
was observed in the parameter region near the continuous 



phase transition 18,22 . In Ref. 22 , by using the gauge fix- 
ing condition, the gauge field propagator is investigated. 
There the gauge field was splitted into a singular (instan- 
ton) part and the remaining regular part. On the criti- 
cal line and also in the confinement phase it was shown 
that propagator of the regular part has a vanishing mass, 
whereas the full propagator has a finite mass. From this 
observation, it is obvious that the finite mass at the crit- 
icality stems from the instanton gas. 

In order to check whether the above observation 18,22 
on the generation of gauge-boson mass on the critical 
line can be applicable also to the present model, one 
needs to separate the effect of nonperturbative instan- 
ton effects. In Ref. 24 the related work was reported for 
the 0(3) sigma model, in which the effect of nonpertur- 
bative spin configuration was studied by suppressing the 
hedgehog configurations. Below we study this problem 
by considering the noncompact U(l) gauge theory of the 
CP 1 Schwinger bosons directly and compare the results 
with those of the compact U(l) system. The action of 
the noncompact model is given by 



5 = -f E (5+^4 + H.c. 



x,fj,,a 



£2 
' 2 



EE( *^) 



(7.13) 



where Q x ,nv is defined by Eq.(6.1) and the constant ei 
controls fluctuations of the noncompact gauge field. For 
small fluctuations of gauge field, |O s , Ml /| <C 1, hence the 
two models are almost equivalent under the relation 62 ~ 

C2- 

In Fig. 28 we present the phase diagram of the model 
(7.13) obtained by calculating the specific heat. The 
crossover line in the compact case disappears as we ex- 
pect from the discussion on the instanton. The region of 
the ordered-Higgs phase is enlarged by the suppression 
of fluctuations of B xll . 




Higgs 



Coulomb 



0.5 1 1.5 2 2.5 3 3.5 4 
FIG. 28. Phase diagram of the noncompact gauge model. 
In the noncompact case, the region of the ordered-Higgs phase 
is enlarged by the suppression of gauge-field fluctuations. The 
crossover line from the dense to dilute instanton regions dis- 
appears. 
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In Fig. 29 we present the gauge-boson mass M' G calcu- 
lated as in the compact case. It is obvious that M' G ~ 
in the disordered phase, and it starts to increase at the 
critical point as expected. Its behavior in the ordered- 
Higgs phase is more or less similar to that of the compact 
gauge model. 

Then we are interested in the instanton density in the 
noncompact model. From the definition of the instan- 
ton density (6.3), it is obvious that instantons can be 
generated even in the noncompact model. In Fig. 31 we 
compare the instanton density in the two models, p{c2) in 
the compact model and p' {e^) in the noncompact model. 
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FIG. 29. Gauge-boson mass in the noncompact gauge 
model as a function of ci. 
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FIG. 31. Instanton density in compact and noncompact 
U(l) gauge theories. The solid curve is a fit p oc exp(— lc?) 



From the results, it is easily seen that p'{c-i) in the non- 
compact gauge model is smaller than pic-i) in the com- 
pact model as expected. However there exists a number 
of instantons for small C2 (i.e., large gauge coupling) and 
one may wonder whether a confinement phase exists for 
sufficiently small e2 . The observation of the specific heat 
and the gauge-boson mass of course denies this possibil- 
ity. 

From the studies on the noncomapct gauge theory in 
this subsection, it is obvious that the observed nonvanish- 
ing gauge-boson mass on the critical line of the comapct 
gauge theory stems from the existence of instantons and 
not the Anderson-Higgs mechanism. 
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FIG. 30. Gauge-boson mass in the noncompact gauge 
model as a function of C2 with ci = 0. 



VIII. DISCUSSION 

In this paper, we have introduced 3D CP 1 +U(1) lattice 
gauge theory and studied its phase structure intensively. 
The specific heat exhibits a second-order phase transi- 
tion at ci = Ci c (c2). The 0(3) spin correlation functions 
show that this critical point separates the spin-ordered 
and disordered phases. The universality class changes 
from the 0(3) spin model at C2 = to the 0(4) spin 
model C2 = 00 smoothly. The MFT interprets the spin 
ordered phase as the Higgs phase and the spin-disordered 
phase as the confinement phase. The distributions of in- 
stantons are studied, which provides us with a rough im- 
age of gauge-field configurations in various points in the 
C2 — ci plane. The gauge-invariant mass Mq of gauge 
boson has two origins; (i) condensation of instantons in 
the confinement phase and (ii) Anderson-Higgs mecha- 
nism (condensation of the CP 1 field) in the Higgs phase. 
As explained above, the conventional MFT plus one-loop 
correction like Eq.(7.7) predicts that Mq = on the crit- 
ical points. However, we observed Mq ^ at ci = c\ c 
due to the remaining instantons. This suggests that the 
system is in the confinement phase just on ci = Ci c (c2). 
Mg(ci c ) is compatible with the second-order phase 
transition because the spin gap Ms{c\ c ) = 0. 

In our recent paper 21 , we introduced a 3D nonlocal 
U(l) lattice gauge theory, which is obtained by mimick- 
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ing the effect of massless (gapless) matter fields in a form 
similar to the gauge model Eq.(7.8). This model con- 
tains an inverse gauge coupling g e , which is propotional 
to the number of massless matter fields. By studying it 
numerically, we found that there exists a deconfincmcnt 
phase transition at a critical gauge coupling g e = g ec . 
This result suggests that a similar deconfinement phase 
may appear on the critical points in the present model 
if the number of matter fields increases as CP 1 — > CP^ 
(iV=large). This possibility is also consistent with the 
recent analytical studies on the massless QED3 which 
suggest that a deconfinement phase is realized for a suf- 
ficiently large number of massless fermions 23 . 

From these motivations, we numerically studied the 
CP N + U(l) gauge theory with N = 2,3 and 4, m par- 
ticular its gauge-boson mass. This model has a qualita- 
tively same phase structure as the CP 1 -|-U(1) model but 
the critical line c\ = c\ c (c<2) shifts upward in the phase 
diagram (shrinkage of the Higgs phase). 



finement phase and the GLT of Eq.(2.17) describes this 
phase transition. 
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FIG. 32. Gauge-boson mass near the critical points in CP^ 
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In Fig. 32, we compare the gauge-boson mass for N — 
1,2,3 and 4. Mq decreases slightly as TV increases but is 
still nonvanishing. We hope to report on more systematic 
study of the CP Ar +U(l) model in a future publication. 

The present study should certainly shed some light on 
rich properties of wide range of quantum spin systems as 
a good Ginzburg-Landau-type phenomenological model. 
For example, the observed phase diagram and the study 
of the instanton in the present CP 1 +U(1) model are quite 
useful to understand the phase structure of the contin- 
uum GLT of Eq.(2.17). As pointed out in Sec.7C, the 
RG flows in Fig. 2 and Fig. 27 are very similar. Further- 
more, we note that the instanton distributions in Fig. 16 
show that the most instantons have topological charges 
\Q X \ less than 3. For the AF Heisenberg model with ring 
exchange, it was argued 14 that instantons with |<3 X | =1,2 
and 3 do not contribute to disordering the gauge dynam- 
ics because of a cancellation mechanism due to the Berry 
phase. Therefore, if there exist the Berry phase term 
in the action of the present model (2.4), the Coulomb 
phase is realized on the critical line instead of the con- 
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